Effect of Crystal-Field Splitting and Inter-Band Hybridization 
on the Metal-Insulator Transitions of Strongly Correlated Systems 
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We investigate a quarter-filled two-band Hubbard model involving a crystal-field splitting, which 
lifts the orbital degeneracy as well as an inter-orbital hopping (inter-band hybridization) . Both terms 
are relevant to the realistic description of correlated materials such as transition-metal oxides. The 
nature of the Mott metal-insulator transition is clarified and is found to depend on the magnitude 
of the crystal-field splitting. At large values of the splitting, a transition from a two-band to a one- 
band metal is first found as the on-site repulsion is increased and is followed by a Mott transition 
for the remaining band, which follows the single-band (Brinkman-Rice) scenario well documented 
previously within dynamical mean-field theory. At small values of the crystal-field splitting, a direct 
transition from a two-band metal to a Mott insulator with partial orbital polarization is found, 
which takes place simultaneously for both orbitals. This transition is characterized by a vanishing 
of the quasiparticle weight for the majority orbital but has a first-order character for the minority 
orbital. It is pointed out that finite-temperature effects may easily turn the metallic regime into a 
bad metal close to the orbital polarization transition in the metallic phase. 

PACS numbers: 71.27.+a,71.70.Ch,71.30.+h,71.10.Fd 



I. INTRODUCTION 

The Mott metal-insulator transition 1 ^ plays a key role 
in the physics of strongly correlated electron materi- 
als. Over the last fifteen years, our theoretical under- 
standing of this phenomenon improved considerably, due 
to the development of the dynamical mean-field theory 
(DMFT) (Refs. i and i). A number of model stud- 
ies were performed in order to clarify the nature of the 
transition in both a single-orbital and multi-orbital con- 

^,{■5, 6, 7,8, 9 , 10.11,12, 13 

In the context of real materials, however, several im- 
portant features must be considered, which are not al- 
ways taken into account in model studies. This includes 
in particular two key aspects: (i) the breaking of orbital 
degeneracy by the crystalline environment and (ii) the 
existence of hopping terms coupling different orbitals on 
different sites of the crystal (inter-orbital hopping or hy- 
bridization). We note at this stage that the breaking 
of orbital degeneracy can correspond to a rather large 
energy scale (of order 1-2 eV) when one has in mind 
the crystal-field splitting between t2 g and e g levels in a 
transition-metal oxide, but it can also correspond to a 
smaller energy scale (a small fraction of an electron volt) 
when considering, e.g., the trigonal splitting of the ti g 
levels induced by a distortion of the cubic symmetry In 
the former case, an effective model with fewer orbitals 
can often be considered, but in the latter case, all or- 
bital components may still be relevant, albeit with dif- 
ferent occupancies, and one has to use a model involving 
several orbitals with slightly different atomic level posi- 
tions. In the present paper, we shall designate the lifting 
of orbital degeneracy by the generic term of "crystal-field 
splitting," but it is mostly the case where this is a small 



energy scale (e.g., trigonal splitting of the tig multiplct) 
that we have in mind for applications. 

Indeed, the physical effects arising from the compe- 
tition of crystal-field splitting and strong correlations 
have attracted a lot of attention recently, in particu- 
lar in LDA+DMFT electronic structure studies of many 
different compounds. We now quote just a few exam- 
ples. Pavarini et al^ pointed out that the lifting of cu- 
bic symmetry by the GdFe03-type distortion plays a key 
role in determining the metallic or insulating characters 
of d 1 transition- metal perovskites such as (Sr/Ca)V03 
(small distortion, metals) and (La/Y)Ti03 (larger distor- 
tion, insulators). Indeed, for a given value of the on-site 
Coulomb repulsion U, the lifting of the orbital degener- 
acy makes the insulating state more easily accessibl e) 9 ' 10 . 
Furthermore, correlation effects considerably enhance the 
effective crystal-field splitting, hence favoring orbital po- 
larization (as also emphasized in Ref. [TH for these com- 
pounds). This correlation-induced enhancement of the 
effective crystal-field splitting and this increased orbital 
polarization have also been show n 16 : 17 to play a key role 
in the metal- insulator transition of V2O3, with the e£ 
component of the ti g level much more occupied than the 
aig component in the insulating phase (see also the pre- 
vious LDA+DMFT studies of V 2 3 in Refs. QJ and 111). 

Such effects were discussed, at a model level, in the pio- 
neering paper of Manini et al&, motivated by the physics 
of fullcrcnc compounds. In this work, a model consisting 
of two orbitals occupied by one electron (quarter-filling) 
was considered, and the combined effect of a crystal-field 
and of on-site repulsion was studied in the framework of 
DMFT. This work identified several phases, most notably 
a two-band metallic phase (with partial orbital polariza- 
tion) and a one-band metallic phase (with full orbital 
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polarization), as well as a fully orbitally-polarized Mott 
insulating phase. 

However, some questions of great importance were left 
unanswered by this early study. To quote just a few of 
these issues: (i) What is the nature of the metal-insulator 
transition in the different ranges of crystal field? (ii) How 
exactly does the crossover between a two-orbital Mott 
transition to a one-orbital Mott transition takes place? 
(Hi) What is the effect of an inter-orbital hopping, al- 
ways present in real materials, and in particular does it 
wipe out the two-band metal to one-band metal tran- 
sition within the metallic phase? and, finally, (iv), is 
it possible to obtain within DMFT the insulating phase 
with partial orbital polarization, which is expected from 
general strong-coupling arguments? (As we shall see, the 
answer is affirmative and this phase was overlooked in 
the DMFT study of Ref. i). 

All these questions are directly relevant to the under- 
standing of real materials (e.g., V2O3 and S^RuC^) and 
to a better qualitative interpretation of the results of 
LDA+DMFT studies. The aim of the present article is 
to provide a detailed answer to these questions. This is 
made possible, in particular, by the recent development 
of numerical techniques for solving efficiently the DMFT 
equations, in particular continuous-time Monte Carlo al- 
gorithmt32ii2^. 

Let us point out that another related model study re- 
cently appeared, namely that of the two-orbital model 
at half-filling (i.e., two electrons in total )£. In this case, 
the physical issues are quite different since one evolves 
from a two-orbital Mott insulator in the absence of crys- 
tal field to a band insulator at large crystal field (not a 
one-orbital Mott insulator as in our quarter-filled case). 
Also, this study did not consider the effect of an inter- 
orbital hopping. In this respect, our work and Ref. [H can 
be considered quite complementary to one another. 

Finally, we emphasize that the interplay between 
crystal-field splitting and strong correlations is made 
even more complex in the presence of Hund's coupling 
and exchange terms. In a study of BaVS3 it was pointed 
out that when Hund's rule wins over crystal-field effects, 
one can observe a compensation between orbital popula- 
tions rather than an enhanced orbital polarizatio n 24 i 25 ' 26 . 
The competition between Hund's coupling and crystal- 
field is also relevant to the physics of cobaltite a 27 ' 28 i 29 ' 30 , 
ruthenates , 31 i 32 ' 33 i 34 and monoxides under pressure 3 -^. In 
the present work, however, we focus on the interplay of 
crystal field and strong correlations and on the nature of 
the Mott transition, in the simplest possible context and 
consider only the effect of an on-site repulsion. 

This paper is organized as follows: In Sec. [Til we m ~ 
troduce the model and some notations. In Sec. IIII A"l we 
present the phase diagram and discuss qualitatively each 
phase. In Sec. IIIIB1 we discuss in details the insulating 
phases, using both an analytical strong-coupling method 
and complete numerical solution of the DMFT equations. 
In Sec. IIII CI we clarify the nature of the various phase 
transitions: from a two-band to a one-band metal and 



from a metal to a Mott insulator, in the different crystal- 
field regimes. Finally, in Sec. IIV[ we consider the effects 
of a finite inter-orbital hopping and also we discuss some 
finite temperature effects in regimes where the two or- 
bitals have very different quasiparticle coherence scales. 



II. MODEL 

We consider a minimal two-band Hubbard model 
with crystal-field splitting and inter-orbital hybridiza- 
tion, given by the Hamiltonian; 
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In these expressions, i is a lattice-site index, k is the mo- 
mentum in reciprocal space, m = 1, 2 is an orbital index, 
and a =|, [ is a spin index. The sum in the interaction 
term runs over all orbital and spin indices except the case 
when m = ml and a — a' and therefore all intraorbital 
and inter orbital Coulomb interactions arc included. A is 
the crystal-field splitting between the two orbitals (A > 
favors the population of the second orbital, m = 2), and 
U is the density-density Coulomb interaction between the 
two orbitals. 

In this article, we focus on quarter-filling (i.e., one elec- 
tron in two orbitals, per lattice site), which is achieved 
by tuning appropriately the chemical potential fi. Wc 
consider only the density-density form of the interaction 
term, and we do not include the Hund's exchange, spin- 
flip, or pair-hopping terms. The motivation for neglect- 
ing these terms is to keep the Hamiltonian as simple as 
possible. Note however that, with one electron per site, 
the effect of these terms is expected to be small and acts 
basically as a renormalization of the on-site U (Refs. [Til 
and [13). 

The kinetic term ifkin is a two-band tight-binding 
Hamiltonian on the three-dimensional cubic lattice (we 
will also use its Bethe lattice counterpart) , which can be 
written (in k space) as 



e(k) V(k) 
V(k) e(k) 



(3) 



where diagonal elements correspond to the simple cubic 
lattice, and the off-diagonal ones have x 2 — y 2 symmetry; 



e(k) 
V(k) 



= 2£(cos k x + cos k y + cos k z ), (4a) 
iky) cos k z . (4b) 



= 2v / 3I^(cos k x — cos k y ) 
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This corresponds to a hopping between identical or- 
bitals on nearest-neighbor sites, equal to —t. The inter- 
orbital hopping connects orbitals m = 1 and m = 2 on 
next-nearest-neighbor sites and is equal in magnitude to 
•v/3V/2. It has a positive sign for the [±1,0, ±1] neigh- 
bors and negative for the [0, ±1,±1] ones. This symme- 
try choice insures that, for all values of V, the on-site 
(k integrated) kinetic Hamiltonian is diagonal in orbital 
space. This is also the case of all local (k integrated) 
quantities in the interacting model, as can be checked by 
expanding the Green's function in power of V. Hence, 
our model is such that the choice of local orbital basis 
set is adapted to the local crystal symmetry. Physically, 
the model [Eq. ([J|] is a reasonable description, for exam- 
ple, of an e g doublet split by the breaking of the cubic 
symmetry. 

For zero hybridization, V=0, the density of states 
(DOS) is reduced to the DOS of the cubic lattice for 
both orbitals shifted by ±A/2. We set the energy unit 
by t — 1/6, or cquivalently D = 1, where D is the half- 
bandwidth. 

We solve this model in the DMFT framework 3 -. Since 
our main aim is to elucidate the nature of the metal- 
insulator transitions in this model, we focus in this article 
on the paramagnetic phases. The self-consistent impu- 
rity problem is solved with two numerical techniques: (i) 
Exact diagonalization (ED) as described in (Refs. d and 
HI), with a "star-geometry" for the bath hybridization 
function using five bath states per orbital degree of free- 
dom; (ii) the recently introduced continuous time quan- 
tum Monte Carlo algorithm (CT-QMC) using an expan- 
sion in the impurity model hybridization functio n 22 ! 23 . 
CT-QMC is more precise than ED and is necessary to 
establish the existence of the partially polarized insulator 
phase (see Sec. IIII 51) . as we shall discuss further below. 



III. RESULTS IN THE ABSENCE OF 
INTER-ORBITAL HYBRIDIZATION 

A. Zero-temperature phase diagram 

The DMFT phase diagram of model [Eq. (P)] at 
quarter-filling and without inter-orbital hybridization 
(V = 0) is presented on Fig. [TJ The effect of a non-zero 
V will be considered in Sec. lIVl The general shape of this 
phase diagram can be easily anticipated by considering 
the various limiting cases^: 

(i) For A = 0, one has a well documented two-band 
degenerate model. The model undergoes a correlation- 
driven Mott transition at a critical U^ =Q ~ 3.76 which 
is close to the results obtained by other authors for 
the Bcthc lattice (semi-circular DOS with identical half- 
width D = 1)&&1. 

(ii) For very large A ^> D, the minority orbital (orbital 
m = 1) is pushed to very high energy and becomes com- 
pletely empty, so that it can be ignored altogether. The 
quarter-filled two-band model thus reduces to a singlc- 




"T" 



T 



FPI 



MIT, V=0 

PPM-FPM line 

Strong coupling, U = 1 / ( 6 A ) 



FPM 



PPM 



0.2 



0.4 0.6 
A AD 



FIG. 1: (Color online) Zero-temperature phase diagram 
(paramagnetic phases), on the cubic lattice without hy- 
bridization (V = 0), and for one electron per site. The solid 
(black) line separates metallic and insulating regions. The 
dot-dashed (red) line separates the partially polarized metal 
(PPM) and the fully polarized metal (FPM) (for details see 
text). The dashed (green) line is the result of the strong- 
coupling mean-field analysis (see Sec. IIII B II and Eq. [TO}: 
It separates the partially-polarized insulator (PPI) from the 
fully-polarized insulator (FPI). Arrows indicate the set of pa- 
rameters used in Figures [4] and [6] The ED solver was used. 



band model at half-filling. This situation has been thor- 
oughly studied within DMFT and yields a correlation- 
induced Mott transition at U^^°° ~ 2.75. (see, e.g., 
Ref. |3j and references therein). The metal-insulator tran- 
sition line (plain/black line on Fig. [TJ) interpolates be- 
tween the limiting critical couplings corresponding to 
A = and A = oo. The system is insulating above 
this line and is metallic below. 

(Hi) The non- interacting model (U = 0) obviously has 
a transition between a two-band metal for A < D and 
a one-band metal for A > D. For A = D, the minority 
band crosses the Fermi level and becomes empty. This 
effective-band transition separating a two-band situation 
at low energy from a non-degenerate band can actually be 
followed through the phase diagram (dashed-dotted/red 
and dashed/green lines on Fig.[T]), as we now discuss. 

We note that we have not attempted to precisely de- 
termine whether the orbital-polarization lines cross the 
metal-insulator transition (MIT) line at a single point, 
or whether the orbital polarization line in the insulating 
phase and in the metallic phase hit the MIT boundary 
at slightly different locations. 

In the absence of hybridization (V — 0), we can use the 
orbital polarization as a faithful indicator of the transi- 
tion between the two-band and a one-band regime. This 
quantity is defined as 



Sn = 



(5) 



\n>) + \n<) 

in which > and < stand for the majority and minority 
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orbitals, respectively, ((n>) > («<))• At quarter filling 
and for A > 0, this reduces simply to Sn = (n 2 — fix). 

As the crystal-field splitting is increased, one reaches 
a critical value at which the orbital polarization reaches 
Sn = 1, indicating a completely empty minority orbital. 
The line along which this happens in the (A, U) plane, is 
indicated by the dashed-dotted (red) line in the metallic 
phase and by the dashed (green) line in the insulating 
phase. Hence, four different phases are apparent on the 
phase diagram of Fig. [TJ a partially polarized (two-band) 
metal (PPM), a fully polarized (one-band) metal (FPM), 
a partially polarized Mott insulator (PPI), and a fully 
polarized Mott insulator (FPI). 

As already pointed out by Manini et aL— , and as clear 
from Fig. [IJ the value of the crystal-field, at which the 
transition from the PPM to the FPM takes place, is 
strongly reduced by interactions. While it is set by the 
half-bandwidth at U = 0, it is renormalizcd down by the 
quasiparticle weight Z> in the presence of interactions. 
Hence, a crystal-field splitting considerably smaller than 
the half-bandwidth can be sufficient to induce a two-band 
to one-band metal transition. 

It is important to realize, however, that the value of A 
needed to fully polarize the system vanishes only in the 
limit U = oo. In other words, the orbitally-degenerate 
Mott insulator at A = has a finite orbital polarizabil- 
ity, even within the DMFT approach. Hence the PPI 
phase at large U and small values of A exists. This point 
was incorrectly appreciated by Manini et al&, largely for 
numerical reasons. Indeed, the ED algorithm is inap- 
propriate to correctly capture the PPI phase. In the 
present article, we establish (Sec. IIII B[) the existence of 
the partially polarized insulating phase within DMFT us- 
ing both an analytical proof at strong-coupling limit and 
a complete numerical solution of the DMFT equations 
based on the new CT-QMC algorithm. 

Let us point out that in this zero-temperature phase di- 
agram, all the transitions are second order, except for the 
transition from the PPM to the PPI - which is second or- 
der for the majority orbital and first order for the minor- 
ity orbital, as will be explained below in Sec. IIII C 21 At 
finite temperatures T > 0, the MIT becomes first order 
throughout the phase diagram, as in canonical DMFT 
solutions, whereas the other transitions remain second 
order. 

In the two following subsections, we describe in more 
details the nature of these different phases and we inves- 
tigate the phase transitions between them. 



B. Existence of the partially polarized insulator 

1. Strong- coupling analysis: Kugel-Khomskii model 

At strong coupling U 3> D (or U 3> t), in the 
Mott insulating phases, an effective low-energy model 
can be derived, following Kugel and Khomskii^l (see also 
Ref. l38l ). The low-energy Hilbert space contains only 



the four states \i, m, a) with one electron on each site 
(m = 1,2; a =\ , J,). The effective Hamiltonian acting on 
these states reads, 

i 

+ {USiSj) + JoiTifj) + J^SiS^if^y (6) 

(«) 

In this expression, (ij) denotes the bonds between 
nearest- neighbor sites, and the spin and pseudo-spin (i.e., 
orbital isospin) operators are given by: 

Si = — ^ ^ dl am T aa idi a i m , (7a) 

m 

Ti = 2 ^ y d ia . m T mm > djerm 1 1 (7b) 



in which f are the Pauli matrices. In particular, the z 
component of these operators (with eigenvalues ±1/2) is 
given by: 

S i = \ ~ fkl2 + fiin ~ n m)> ( 8a ) 
T* = \{hi\2 + «U2 - hi\\ - fni\). (8b) 

The (superexchange) couplings J s , J , and J m are given 
by^ 



J 2t 2 

Js = Jo = — = — = J- 



(9) 



The particular symmetry between these couplings is due 
to the choice of a density-density interaction and to the 
neglect of the Hund's exchange. 

At strong coupling, in the insulating phase, the DMFT 
solution of the original model [Eq. (JT])] reduces to a 
static mean- field solution of Eq. ©. Focusing on the 
non-magnetic phase ((Sf) = 0), the orbital polarization 
Sn = 2(T Z ) is given by the self-consistent equation, at 
finite temperature T = 1//3: 



where 



Sn = tanh ^— (A — A c 5n)^, 



zJ t 2 
A c = — = z—, 
2 U 



(10) 



(11) 



is a critical value of the crystal-field splitting and z 

is the coordination number of the lattice (number of 

nearest neighbors). For the simple cubic lattice, with 

z = 6 and half-bandwidth D = zt, this yields, A^ L = 

D 2 /(6U), while for the large-connectivity Bethe lattice 

with nearest-neighbor hopping t = D/(2y/z), one has: 
a bl = D 2/(4uy 
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At zero temperature (f3 = oo), the solution of Eq. ([XDI 
reads, 



Sn 




A < A c 
A > A r 



(12) 



Hence, this shows that the orbitally degenerate insulator 
has a finite orbital susceptibility at T = 0, Xorb = 1/A C , 
and that a finite crystal-field A = A c must be applied to 
fully polarize the insulating phase. The strong-coupling 
expression A = A c = D 2 /(6U) for the cubic lattice cor- 
responds to the dashed (green) line displayed on Fig. [U 
separating the PPI from the FPI phases at T = 0. 

At finite temperature, a good approximation to the 
solution of Eq. ([TO)) turns out to be; 



Sn 



A 



1 



(13) 



Finally, we would like to emphasize that, when think- 
ing of DMFT as an exact method in the limit of large 
lattice coordination z — ► oo, it is quite clear that a non- 
zero value of the critical A c (and hence a finite extent of 
the PPI phase) is to be expected. Indeed, the orbital ex- 
change coupling [Eq. (JO])] scales as 1/z (since t oc l/^fz), 
hence, the critical A c is of the order of the exchange field 
between a site and all its neighbors, i.e., of order zJ, 
which remains 0(1) as z — * oo. The uniform orbital sus- 
ceptibility of the orbitally degenerate Mott insulator is 
indeed finite at T = (this should not be confused with 
the fact that the local susceptibility would scale as 1/z 
and hence vanish in the large- z limit). 



2. Numerical solution: Importance of global moves in the 
quantum Monte Carlo algorithm 

The analytical estimate at finite temperature [Eq. (p~5]) ] 
provides a very useful benchmark when solving numeri- 
cally the DMFT equations for the original model in the 
strong-coupling regime. Indeed, it is actually non-trivial, 
from the numerical point of view, to successfully stabi- 
lize the partially polarized insulating phase. To achieve 
this, we have used the CT-QMC method, and it proved 
necessary to implement global Monte Carlo moves, in 
addition to the Monte Carlo moves proposed in Ref. 8. 
In CT-QMC, a configuration is given by a collection of 



fermionic operators c ai (nj 



,(tjv) at different imag- 



inary times r, and the a, are the fermionic species of the 
operators. The global moves arc implemented by chang- 
ing all (Xi into a new set of o! i and accepting the move 
with a probability satisfying the detailed balance condi- 
tion. In this work, we have used two global moves that 
switch the spin (T^l) and the orbital (1 <-> 2) indices. In 
the absence of these global moves, the calculation can be 
trapped in some regions of the phase space at low tem- 
perature, leading to a wrong (overestimated) value of the 
polarization. 




CT-QMC wgm 
CT-QMC w/o gm 
Strong coupling 



FIG. 2: (Color online) Temperature dependence of the orbital 
polarization, 8n, for the cubic lattice in the PPI phase. Black 
dots are the CT-QMC data with the use of global moves (gm) 
in the spin and orbital space (see text for details). (Cyan) 
squares are the CT-QMC data without global moves. The 
solid (green) line is the strong coupling result given by Eq. 1101 
The arrow shows the zero-temperature value of the polariza- 
tion in the strong coupling limit, Sn = 0.84. These results are 
obtained for U=4, A=0.035 and using the CT-QMC solver. 



This is illustrated in Fig. [2l which displays the tem- 
perature dependence of the polarization in the insulating 
phase, at small A. The result of Eq. (fTUj) is compared 
to the CT-QMC results with and without global moves. 
One can see that without global moves, the polarization 
is bigger than its strong-coupling value, whereas the con- 
trary is expected. This gives a clear indication that global 
moves are needed. When the correct implementation of 
the CT-QMC algorithm with global moves is used, the 
polarization falls below its strong-coupling value. Note 
that these results are actually obtained for an interme- 
diate value of U = 4, which shows that the range of 
validity of the strong-coupling approximation is actually 
quite extended. The agreement between the DMFT data 
with global moves and the strong-coupling result is seen 
to be excellent and both indubitably show the existence 
of the partially polarized insulator. 

We have not been able (as in Ref. @), when using the 
ED solver at T = in the insulating phase, to stabi- 
lize the partially polarized insulating solution at small A. 
This is probably because this solution is too delicate and 
involves a number of competing low-energy scales (J, A) 
to be faithfully reproduced given the simple parametriza- 
tion and limited number of states in the effective bath, 
which can be handled within ED in a two-orbital context. 
However, ED performs quite well in the metallic phase, 
and it is quite instructive to compare the iso-polarization 
lines (Sn — const.) in the (A, U) plane, determined from 
ED, close to the metal-insulator transition, to the strong- 
coupling result A/A c = Sn (i.e., UA = zt 2 8n). This 
comparison is made in Fig. [3] in the case of the Bcthe 
lattice (for simplicity) . The ED data on the metallic side 
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FIG. 3: (Color online) Metal-insulator transition (red) and 
iso-polarization lines for the Bethe lattice. Different colors 
mark different values of the polarization (see legend). Dashed 
lines are ED results, while solid lines are the solutions of 
Eq. [12] for constant polarization. 

of the transition match very well to the strong-coupling 
form of the iso-polarization lines on the insulating side. 
Thus, this provides a complementary way, starting from 
the metal, to document the existence of the PPI regime. 



C. Metallic phases and the nature of the 
metal-insulator transition 

We now turn to the metallic phases. There, the self- 
energies can be Taylor expanded at low-frequency as 

RE^(w + i0 + ) = SRS^(0) + (1 - 1/Z^)u + ■■■ , (14) 

in which Z > and Z < are the quasiparticle weights of the 
majority and minority bands, respectively. The quasi- 
particle weights, Z^s = (1 — d^'E^(iui)/diu)) l^^o were 
extracted from the imaginary frequency data with the 
use of third-order polynomials. The minority and major- 
ity Fermi surfaces in the metallic phase PPM are deter- 
mined, respectively, (for V = 0) by: 

e(k) = M +^-5RE > (0) = ,u > , (15a) 

e(k) = M -^-5RE<(0) = M< . (15b) 

The quantities /i>,/i< can be viewed as effective crystal- 
field levels renormalized by interactions (or effective 
chemical potentials for each type of orbitals) , and a renor- 
malized crystal-field splitting can also be defined as 

A off = A + 5RS<(0)-5RS>(0). (16) 

The various transitions arc conveniently described in 
terms of and Z>. On general ground, there are two 



simple mechanisms by which a given orbital can undergo 
a transition from a metallic behavior to an insulating one: 

(i) The quasiparticle weight Z may vanish at the MIT. 
This is the well-known Brinkman-Rice scenario, which 
is realized, e.g., within the half-filled single-band Mott 
transition within DMFT. It is also realized for degenerate 
orbitals with A = 0: Z > = Z < vanishes continuously at 
U C A =°- 

(ii) It may also happen that cither of the equations 
[Eq. (fT5|) ] fails to yield a solution, i.e., the "effective chem- 
ical potentials" /i> or /i< move out of the energy range 
[— D, +D] spanned by e(k). This, in turn, can happen in 
a continuous or in a discontinuous way. 



1. Orbital polarization and metal-insulator transitions at 
large crystal field 

We first consider values of the crystal-field splitting 
larger than > 0.1. Two successive transitions are ob- 
served as U is increased, from a two-band metal (PPM) 
to a single-band metal (FPM) - followed by a metal- 
insulator transition (FPM to FPI) . Figure H] (top panel) 
displays the quasi-particle residues Z > , Z < and orbital 
polarization 5n as U is increased at a fixed A = 0.3 (in- 
dicated by the arrow on Fig. [T]) . The lower panel of Fig. 0] 
displays ii^ and A c g. 

For U < C/p(~ 2.2), in the two-band metallic phase 
(PPM), both quasiparticle weights decrease as U is in- 
creased, and the orbital polarization gradually increases. 

At U = Up, the polarization saturates to 5n = 1 and 
the minority band becomes empty. This happens follow- 
ing the mechanism (ii) above: the minority band effective 
level position hits the bottom of the band (/x< = — D 
at U = Up) and the renormalized crystal-field splitting 
reaches A e g = +D [as clearly seen from Fig. 2] (lower 
panel)]. Simultaneously, /i> vanishes at Up and remains 
zero for U > Up. This indicates that particle-hole sym- 
metry is restored at low-energy for the majority band 
throughout the FPM phase. 

For Up < U < Umit, the minority band is empty and 
becomes inactive. The remaining half-filled majority or- 
bital forms a single-band metal and is subject to the local 
Coulomb interaction. This is illustrated in Fig. [5] where 
we plot the DOS of both orbitals. Note that the major- 
ity orbital very quickly becomes particle-hole symmetric 
over its full bandwidth as U increases. The quasi-particle 
weight of the majority band Z > is strongly reduced in 
this regime. Note that neither Z > nor Z < vanishes at 
the orbital polarization transition Up. In fact, also the 
minority (empty) band sclf-cncrgy remains linear in fre- 
quency at low energy in this regime, and a Z < can still 
be formally defined (as plotted on Fig. 2]), although it no 
longer has the physical meaning of a quasiparticle spec- 
tral weight since there is no Fermi surface for that band. 
In particular, the increase of Z < in this region should not 
be interpreted as a decrease in the correlation effects. 

Eventually, the transition from a single-band strongly 
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FIG. 4: (Color online) PPM-FPM-FPI transitions along the 
constant A = 0.3 line for the cubic lattice without hybridiza- 
tion (V = 0). 

Top panel: Orbital polarization, 5n and QP residues, Z > and 
Z< are shown by (black) dots, filled (green) diamonds and 
open (red) squares, respectively. 

Bottom panel: Effective crystal field splitting, A e // and ef- 
fective chemical potential for both bands, ^i > , fi < are repre- 
sented by (black) dots, filled (green) diamonds and open (red) 
squares, respectively. 

The vertical lines show the full polarization (violet) and MIT 
transitions (magenta). The horizontal (brown) lines show the 
top and bottom of the bare band. The ED solver was used. 



correlated metal to a Mott insulator with full orbital po- 
larization is found at U = U U n t/ c A= °° ^ 2.75). 
The nature of this transition has been exhaustively de- 
scribed in the context of DMFT studies of the single-band 
model: Z > vanishes continuously at the critical point and 
the metal- insulator transition is second order (at T=0). 
The low- frequency majority self-energy 9?£>(w + i0 + ) 
acquires a pole on the real frequency axis in the insu- 
lating phase. The location of this pole depends on the 
choice of the chemical potential within the insulating gap. 
For a specific choice (as done in Fig. [3]), the pole is lo- 
cated at zero-frequency so that particle-hole symmetry 
is restored at low energy and the self-energy diverges as 
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FIG. 5: (Color online) Density of states in the 1-band metallic 
phase. The majority (minority) orbital is shown in the upper 
(lower) half of the plot. The ED solver was used with U = 2.4 
and A=0.3. 

9f?£>(cj + i0 + ) ~ l/cu. 

It should be emphasized that the very small value of 
Z > in the one-band (FPM) metallic phase implies that 
the quasiparticles are actually quite fragile in that phase 
and can be easily destroyed by thermal effects. Hence, 
the orbital polarization transition at U = Up from a 
two-band to a one-band metal at T = may actually 
appear, at finite-temperature, as a transition between a 
two-band metal and a one-band incoherent "bad metal" 
(or quasi- insulator). We shall come back to this point in 
more details in Sec. ITVBl 



2. Metal-insulator transition at small crystal field 

In the small crystal-field regime (A < 0.1), to the best 
of our numerical accuracy, there appears to be a simul- 
taneous metal-insulator transition for both orbitals from 
a two-band metal (PPM) to a Mott insulator with par- 
tial orbital polarization (PPI). Note that in this region 
we needed to have recourse to finite-temperature Monte 
Carlo simulations. 

The nature of the MIT has been well documented by 
previous DMFT studies in the degenerate case A = 0. At 
T = 0, the transition is second order with a quasiparticle 
weight Z > = Z < vanishing continuously at U^ =0 , while 
at T > this transition is first order. 

On the top panel of Fig. [51 we display the quasipar- 
ticle weights Z >1 Z < as a function of U for /3 = 100 
and for a small value of A = 0.05, along with the or- 
bital polarization Sn. The MIT takes place at a criti- 
cal coupling U^, which is smaller than U^ =0 (Fig. [I}. 
Note that the data in Fig. [5] is obtained for a finite tem- 
perature and, therefore, the critical is also smaller 
than its zero-temperature counterpart (shown in Fig. [I}. 
The orbital polarization continuously increases with the 
interaction and does not approach the value 5n—l at 



8 



SJ A 0.8 h 

. v 

N 0.6 - 



0.4 - 



2 - 



-2 -. 



• — • Sn 

— Stroilg coupling 





, A = A - ReZ + ReL 

eff > < 

i (l = \i. - A/2 - Re"L < 
u = (i + A/2 - ReE > L 



i ■ i ■ i ■ i 



1 1.2 1.4 I.I 
(l/D 



1 



2 3 
U/D 



After the transition, the chemical potential /i can be 
placed (at T = 0) anywhere within the charge gap, and 
therefore, the effective chemical potentials [Eq. (fT5|) ] are 
not longer defined in a unique manner. As documented in 
previous work^i on the orbitally degenerate case within 
DMFT, we expect the majority orbital self-energy to have 
a pole on the real frequency axis, at a position that de- 
pends on /j. For a special choice of fi, this pole is lo- 
cated at uj = 0, which should correspond to a divergence 
of £>(w = 0) and to a divergent self-energy ~ 1/ui at 
low frequency. In order to document this behavior, we 
plot in the inset of Fig. [5] the real part of the imaginary 
frequency self-energies, 3?S^(i7r//?) at the first Matsub- 
ara point as a function of n (for a given value of the 
interaction U = 3). One can clearly see that the major- 
ity orbital self-energy, 5RE>(i7r//?) becomes very big and 
changes sign at fj, ~ 1.22 while 5RS < (i7r//3) stays con- 
stant within the gap. A careful scaling analysis shows 
that S>(w = 0) indeed diverges at a critical value of 
fi. Together with the vanishing of Z>, this shows that 
the transition for the majority orbital follows the stan- 
dard DMFT scenario identical to the orbitally degener- 
ate case. Furthermore, since 5RS<(i7r//3) does not vary 
significantly when /i is varied within the gap, one can 
unambiguously define also in the insulating phase. In 
contrast, depends on the choice of fj.. One should 
note here that the chemical potential, (i, defined in this 
way in the insulating phase continuously connects to the 
chemical potential in the metallic phase. 



FIG. 6: (Color online) PPM-PPI transition along the con- 
stant A=0.05 line for the cubic lattice without hybridization. 
Top panel: The (black) solid dots show the orbital polariza- 
tion, 8n. The QP residues, Z > and Z < , are shown by the 
(green) filled diamonds and (red) open squares, respectively. 
The (blue) dashed line is the strong coupling result. 
Bottom panel: Effective crystal field splitting, A e // and ef- 
fective chemical potential for both bands, ^i > , fi < are repre- 
sented by (black) dots, filled (green) diamonds and open (red) 
squares, respectively. The inset shows the real part of the self- 
energies at the first Matsubara frequency, (in / f3), versus 
chemical potential within the gap for (7 = 3. 
The vertical magenta line shows the MIT. The CT-QMC 
solver was used with (3 = 100. 

the transition point. The minority orbital quasiparticle 
weight Z < remains larger than the majority one Z > in 
the metallic phase. Although it is a delicate issue numer- 
ically, our data appear to be consistent with a majority 
orbital quasiparticle weight Z > , which vanishes continu- 
ously while Z< remains finite at the transition. Note that 
both the majority and minority orbital effective chemi- 
cal potentials [Eq. (|15[) ] stay well within the energy band 
[-D, +D] for all couplings in the metallic phase. The 
transition into the insulating phase for the minority or- 
bital takes place by having n < jumping out of the energy 
band in an apparently discontinuous manner, as we now 
describe in more details. 



In Fig. [6] (bottom panel) , we display these two quan- 
tities as a function of U, choosing for the special value 
at which £> behaves as 1/lu at low frequency. From this 
plot, we see that the minority band becomes insulating 
because /i< is jumping out of the energy band in a man- 
ner that appears as discontinuous (up to our numerical 
precision). Hence, in contrast to the orbital polarization 
transition of the large A case described above, the MIT 
at small A appears to occur in a discontinuous manner, 
as far as the minority band is concerned, while being 
continuous (Brinkman-Rice like) for the majority band. 
Note also that the minority-orbital self-energy has a lin- 
ear behavior at low frequency throughout the insulating 
phase. 



Note that in this finite-temperature calculation, the 
orbital polarization never reaches Sn = 1 as U is further 
increased. From the strong-coupling calculation, we ex- 
pect that it will saturate at Sn ~ 0.987 when U — > oo. At 
zero temperature, however, there is a second-order tran- 
sition at a finite critical value of U where the polarization 
reaches Sn — 1. 
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FIG. 7: (Color online) Zero-temperature phase diagram of 
the cubic lattice with hybridization V = 0.07 and for one 
electron per site. The (black) solid line separates metallic 
and insulating regions. The (black) dot-dashed line separates 
the 2-band and 1-band metals. For the sake of comparison, 
the corresponding zero-hybridization (V = 0) lines are shown 
(in green). The ED solver was used. 



t 1 1 1 1 r 




FIG. 8: (Color online) Orbital polarization 5n (dots/black), 
quasiparticle weights Z2 (green/filled diamonds) and Z\ 
(red/open squares), as a function of U for a fixed value of 
A = 0.3 and a finite inter-orbital hybridization V = 0.07. For 
the sake of comparison, the orbital polarization for V = 
is also displayed (cyan/dashed line). The vertical (magenta) 
line shows the MIT. The ED solver was used. 



IV. EFFECT OF AN INTER-ORBITAL 
HYBRIDIZATION V(k) 

A. Low-energy effective-band transition 

In this section, we consider the effect of a finite hy- 
bridization (inter-orbital hopping V^(k) ^ 0). At low 
values of A, the metal-insulator transition is pushed to 
higher values of U when turning on a small V. While at 
larger values of A, the MIT line is less sensitive to V (as 
illustrated on Fig. [7]). This is expected since at low A the 
inter-orbital hopping increases the kinetic energy in both 
bands while at higher A the hybridization with a band, 
which is already empty, has a smaller effect on the crit- 
ical coupling. As we will discuss in more details below, 
in the presence of the hybridization, the fully polarized 
phases (FPM and FPI) disappear. However, there is still 
a transition from a two-band to a one-band metal at low 
energy. This transition line is pushed up at low values 
of the crystal-field splitting because of the increase in ki- 
netic energy. In non- interacting limit, the finite value of 
V acts as a k-dependent enhancement of the crystal field 
A, and therefore, at small values of the interaction, the 
two-band to one-band transition line is below the corre- 
sponding V = line. 

One should note that the majority (minority) band 
does not have a unique two (one) orbital character, and 
the band index > (<) has to be distinguished from the 
orbital index two (one). 

On Fig. [51 we display the quasiparticle weights and 
orbital polarization as a function of U, for a fixed value 
of V and a rather large crystal field A = 0.3. One clearly 
sees that the MIT follows a similar mechanism than in 
the V = case: only Z 2 vanishes continuously at the 



transition, while Z\ is always finite. 

A noticeable difference with the V = case is that 
the orbital polarization 5n = n 2 — ri\ does not reach 
saturation (Sn < 1) before the MIT (Fig. [8j. This is 
expected, because the low-energy bands in the metallic 
state no longer have a unique (1, 2) orbital character, as 
we now discuss. 

In order to understand more precisely the nature of 
the metallic phase, we use the low-frequency expansion 
of the self-energies and we obtain the expressions of the 
low-energy majority and minority bands, which read, 

2w<(k) = (17a) 

= Ziei k + Z 2 e 2k + (Ziei k - Z 2 e 2 k) 2 + 4ZiZ 2 V k 2 , 
2w>(k) = (17b) 

= Z 1 e lk + Z 2 e 2k - yj (Zie lk - Z 2 e 2k ) 2 + AZ 1 Z 2 V^. 

In these expressions eik = ey — /1 + A/2 + 5RSi(0) and 
£2k = e k - n - A/2 + 5RS 2 (0). The Fermi surface (set 
by lo = 0) is determined by the following condition (in 
which the weights Z\ %2 do not appear): 

= eike2k-V k 2 (18) 
= [e k -M + A/2 + »Ei(0)] x 
[e k -M- A/2 + 9fffi 2 (0)]- Vk 2 . 

We recall that, when V = 0, an orbital polarization tran- 
sition is first encountered at U = Up , at which the Fermi- 
surface sheet corresponding to orbital 1 (determined by 
eik = 0) disappears, since fi — A/2 — 5RSi(0) reaches 
the band-edge. In the presence of V 7^ 0, a similar phe- 
nomenon occurs for the minority low-energy band uj < (k) : 
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FIG. 9: (Color online) Linearized band structure along symmetry lines of the cubic lattice for different values of the interaction. 
We used A=0.3, V—0.07 and the ED solver. Fatness shows a contribution of the spectral weight of the less occupied orbital 
(m = 1) to the majority band (see text for details). [7=2.4 corresponds to the value where the effective crystal field splitting 
exceeds the bare bandwidth and the physical picture effectively becomes single band. The rightmost panel corresponds to the 
PPI solution and we used 5RE2(w + i0 + ) = 5RE2(0) + 0,/oj for the divergent orbital. 



one of the two sheets, which constitute the solution of 
Eq. (fl"8)) ceases to exist. This is expected from conti- 
nuity arguments in view of Eq. (|18|) and of the situa- 
tion at V = 0. This is furthermore demonstrated by 
Fig. [5J which displays the majority and minority low- 
energy bands tu^ (k) along the main directions in the Bril- 
louin zone, as U is increased. It is clearly seen from this 
figure that for U ~ 2.4 (before the MIT, which takes place 
at U ~ 2.8), an effective band transition occurs between 
the two-band metal and a one-band metal at low energy. 
The critical coupling for this effective-band transition is 
slightly increased as compared to its value at V = 0. 

For V = 0, the majority eigenstate |k, >) (correspond- 
ing to eigenvalue u>>(k)) has a unique orbital character 
to = 2. In contrast, for V ^ 0, it has a component on 
both orbital 2 and orbital 1. As a result, the orbital 
polarization does not reach Sn = 1 (Fig. |SJ) at the effec- 
tive band transition between a two-band and a one-band 
metal. On Fig. [5J we have used a "fat band" represen- 
tation to illustrate this point: at each k point, we plot a 
bar whose extension is proportional to the matrix clement 
|(l|k, >)| 2 , measuring the projection of the less-occupied 
orbital to = 1 onto the majority band. 

As U is increased beyond the effective-band transition, 
one is left with a single effective low-energy band, char- 
acterized by a quasiparticle weight, 

z>w = [ei *i e ? )z f , d9) 

where k lies on the Fermi surface of the majority band [sc 
Eq. (US))]. The subsequent Mott transition is character- 



ized by a vanishing quasiparticle weight for the majority 
band Z > ~ Z2 — > 0, as clearly seen from Fig. [8] and from 
the narrowing of that band in the third panel of Fig. [5J 
The key conclusion of this section is that, even in the 
presence of a finite inter-orbital hopping, two distinct 
transitions are observed as U is increased (in the large 
crystal-field regime): first, a second-order transition from 
a metal with two active bands at low energy to a metal 
with only one active band at low energy, and followed by 
a Mott metal-insulator transition of the one-band type. 



B. Orbital-selective coherence and the two-band 
metal to one-band bad-metal transition 

We have seen above that, in a rather extended region 
of the metallic phase, the quasiparticle weight of the ma- 
jority orbital is much smaller than that of the minority 
one. This is especially true close to the two-band to one- 
band metal transition, where Z2 <C Z\. This implies 
that thermal effects can easily destroy the fragile quasi- 
particles of the majority band. This has physical con- 
sequences, which may be important in practice. For ex- 
ample, the two-band metal to one-band metal transition 
at finite temperature may appear in practice as a quasi- 
mctal-insulator transition or more precisely as a transi- 
tion between a two-band metal and a bad (or incoherent) 
metal. This will happen when the temperature, at which 
the system is studied, is higher than the (small) quasi- 
particle coherence temperature of the majority band. 

In order to illustrate this point, we performed finite- 
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FIG. 10: (Color online) Top panel: Imaginary part of the self- 
energies, 3Ei(itj n ) and 3E2(iw n ) for different temperatures 
(see legend for temperature coding). The inset shows the ex- 
trapolation to zero of the imaginary part of the self-energies, 
Q ! Ei i 2(0) versus temperature. 

Bottom panel: Imaginary part of the Green's functions, 
QGi(iu> n ) and SsG2(iw„) for different temperatures (the color- 
coding is the same). The inset shows the density of states at 
the chemical potential pi, 2(0) versus temperature. We used 
U = 2, A = 0.3, V = 0.07 and the CT-QMC solver. 



temperature studies for the following parameter values: 
U = 2, A = 0.3, and V = 0.07, which correspond to 
the two-band metallic regime, not very far from the two- 
band to one-band metal transition. For these parameters, 
the two quasi particle residues arc Z2=0. 34 and Zi=0.59 
(see Fig. |8J). In Fig. [TTJ we display the imaginary part 
of the Green's functions (bottom) SGi^f^Wn) and self- 
energies (top) QT,i^(itOn) on the Matsubara axis, for dif- 
ferent temperatures. In the insets of this figure, we dis- 
play the extrapolated zero-frequency value 3Si i 2(i0 + ), 
which is related to inverse quasiparticlc lifetime and zero- 
frequency density of states, p 1.2(0) = — 3Gi.2(i0 + )/7r, 
respectively. 

It is seen from these figures that, while the minority or- 
bital quantities have quite little temperature dependence, 
the majority orbital, in contrast, displays very strong 



temperature dependence. For example for T > 0.03 (i.e., 
a rather low-energy scale as compared to the bandwidth) , 
the majority orbital is clearly incoherent with a small 
quasiparticlc lifetime and much reduced value of the local 
density of states. At those temperatures, the frequency 
dependence of the self-energy is clearly non-metallic, ex- 
trapolating to a large value at zero frequency. Only at 
a low temperature T ~ 0.01 (200 times smaller than the 
bandwidth), the behavior of a coherent metal is recov- 
ered, with a linear Matsubara frequency dependence of 
^sTi2(iu> n ) extrapolating to a small value at low frequency 
(corresponding to a large quasiparticlc lifetime). 



V. CONCLUSION 

In this paper, we have investigated how a crystal-field 
splitting, by lifting orbital degeneracy, affects the Mott 
metal-insulator transition in the presence of strong on- 
site correlations. The study was performed on a sim- 
ple two-orbital model at quarter filling (one electron per 
site), and we have also considered the effect of an inter- 
orbital hopping (hybridization), which is important for 
applications to real materials. 

Within the metallic phase, a second-order transition 
from a two-band to a one-band metal takes place as the 
crystal field is increased. The critical value of the crystal- 
field splitting, at which this transition takes place, is con- 
siderably lowered for strong on-site repulsion (i.e., the 
effective crystal- field splitting is considerably enhanced). 
This transition has the nature of a effective band tran- 
sition for the rcnormalizcd low-energy bands (i.e., the 
minority band is pushed up in energy and does not cross 
the Fermi energy anymore) and survives in the presence 
of an inter-orbital hopping. 

The nature of the Mott metal-insulator transitions in- 
duced by on-site repulsion was found to depend on the 
magnitude of the crystal-field splitting. At high enough 
values of this splitting, the Mott transition is between 
a one-band metal and a one-band Mott insulator (con- 
ventional Brinkman-Rice scenario): only the majority 
orbital is involved, and the transition is second order 
and characterized by a vanishing quasi-particle weight 
for that orbital. At low values of the crystal-field split- 
ting, the transition is from a two-band metal to a Mott 
insulator with partial orbital polarization. It takes place 
simultaneously for both orbitals: although the transi- 
tion is still continuous for the majority orbital, it has 
a first-order character for the minority orbital. Elucidat- 
ing these transitions and, in particular, establishing the 
existence of the partially orbitally polarized Mott insula- 
tor at low crystal fields was made possible by the recent 
development of the CT-QMC algorithm for the solution 
of the DMFT equations. 

If a finite hybridization (V 7^ 0) is taken into ac- 
count, it is no longer possible to fully polarize the sys- 
tem. Therefore, the FPM and the FPI phases disappear. 
However, there is still a transition from a two-band to a 
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one-band metal at low energy so that the introduction 
of a finite V does not modify the overall picture of the 
model. 

We have also studied the influence of the tempera- 
ture on the two-band metal just below the transition to 
the one-band metal. The temperature can easily drive 
the system into a regime where the quasiparticle weight 
of the majority band is destroyed and the system effec- 
tively becomes a single-band metal. Further increase of 
the temperature above the characteristic temperature of 
both bands leads the system into an incoherent (or bad) 
metal. 

Our study has direct relevance for the interpretation of 
the metal-insulator transitions of transition-metal oxides 
(see Sec. HJ, often accompanied by an enhanced orbital 
polarization. 



Acknowledgments 



We are very grateful to O.K. Andersen, S. Bier- 
mann, A. Rubtsov, and A. Lichtenstein for the discus- 
sions related to this work. We also thank V. Anisi- 
mov, F. Lcchermann, A. Millis, and P. Werner for the 
useful conversations. We acknowledge the support from 
CNRS, Ecole Poly technique, the Agcncc Nationale de la 
Recherche (under contract ETSF), and the Marie Curie 
Grant No. MIF1-CT-2006-021820. This work was sup- 
ported by a supercomputing grant at IDRIS Orsay under 
Project No. 071393 (for the ED results) and at CEA- 
CCRT under Project No. p588 (for the CT-QMC calcu- 
lations). 



1 N. Mott, Metal-Insulator Transitions (Taylor and Francis 
Ltd., London/Philadelphia, 1990). 

2 M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

3 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

4 G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. Parcollet, and C. A. Marianetti, Reviews of Modern 
Physics 78, 865 (2006). 

5 T. Pruschke and R. Bulla, The European Physical Journal 
B 44, 217 (2005). 

6 N. Manini, G. E. Santoro, A. Dal Corso, and E. Tosatti, 
Phys. Rev. B 66, 115107 (2002). 

7 M. J. Rozenberg, Phys. Rev. B 55, R4855 (1997). 

8 P. Werner and A. J. Millis, Phys. Rev. Lett. 99, 126405 
(2007). 

9 E. Koch, O. Gunnarsson, and R. M. Martin, Phys. Rev. B 
60, 15714 (1999). 

S. Florens, A. Georges, G. Kotliar, and O. Parcollet, Phys. 
Rev L B 66, 205102 (2002). 

11 Y. Ono, M. Potthoff, and R. Bulla, Phys. Rev. B 67, 
035119 (2003). 

12 J. E. Han, M. Jarrell, and D. L. Cox, Phys. Rev. B 58, 
R4199 (1998). 

13 L. Laloux, A. Georges, and W. Krauth, Phys. Rev. B 50, 
3092 (1994). 

14 E. Pavarini, S. Biermann, A. Poteryaev, A. I. Lichtenstein, 
A. Georges, and O. K. Andersen, Phys. Rev. Lett. 92, 
176403 (2004). 

15 M. Mochizuki and M. Imada, Phys. Rev. Lett. 91, 167203 
(2003). 

16 M. S. Laad, L. Craco, and E. Miiller-Hartmann, Phys. Rev. 
Lett. 91, 156402 (2003). 

17 A. I. Poteryaev, J. M. Tomczak, S. Biermann, A. Georges, 
A. I. Lichtenstein, A. N. Rubtsov, T. Saha-Dasgupta, and 
O. K. Andersen, Phys. Rev. B 76, 085127 (2007). 

18 G. Keller, K. Held, V. Eyert, D. Vollhardt, and V. I. Anisi- 
mov, Phys. Rev. B 70, 205116 (2004). 

19 M. S. Laad, L. Craco, and E. Muller-Hartmann, Phys. Rev. 



B 73, 045109 (2006). 

20 A. Rubtsov and A. Lichtenstein, JETP Letters 80, 61 
(2004). 

21 A. N. Rubtsov, V. V. Savkin, and A. I. Lichtenstein, Phys. 
Rev. B 72, 035122 (2005). 

22 P. Werner, A. Comanac, L. de' Medici, M. Troyer, and 
A. J. Millis, Phys. Rev. Lett. 97, 076405 (2006). 

23 P. Werner and A. J. Millis, Phys. Rev. B 74, 155107 (2006). 

24 F. Lechermann, S. Biermann, and A. Georges, Phys. Rev. 
Lett. 94, 166402 (2005). 

25 F. Lechermann, S. Biermann, and A. Georges, Progress Of 
Theoretical Physics Supplement 160, 233 (2005). 

26 F. Lcchermann, S. Biermann, and A. Georges, Phys. Rev. 
B 76, 085101 (2007). 

27 C. A. Marianetti, G. Kotliar, and G. Ceder, Phys. Rev. 
Lett. 92, 196405 (2004). 

28 H. Ishida, M. D. Johannes, and A. Liebsch, Phys. Rev. 
Lett. 94, 196401 (2005). 

29 C. A. Perroni, H. Ishida, and A. Liebsch, Phys. Rev. B 75, 
045125 (2007). 

30 C. A. Marianetti, K. Haule, and O. Parcollet, Phys. Rev. 
Lett. 99, 246404 (2007). 

31 A. Liebsch and A. Lichtenstein, Phys. Rev. Lett. 84, 1591 
(2000). 

32 V. I. Anisimov, I. A. Nekrasov, D. E. Kondakov, T. M. 
Rice, and M. Sigrist, European Physical Journal B 25, 
191 (2002). 

33 A. Liebsch and H. Ishida, Phys. Rev. Lett. 98, 216403 
(2007). 

34 X. Dai, G. Kotliar, and Z. Fang (20 06), URL 
http : //arxiv . org/abs/cond-mat/061107 5vl| 

35 J. Kunes, A. V. Lukoyanov, V. I. Anisimov, R. T. Scalet- 
tar, and W. E. Pickett, Nature Materials 7, 198 (2008). 

36 M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 

(1994) . 

37 K. Kugel and D. Khomskii, Sov. Phys. Usp. 25, 231 (1982). 

38 D. P. Arovas and A. Auerbach, Phys. Rev. B 52, 10114 

(1995) . 



